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1.  Overview 

This  work  has  resulted  in  a  number  of  publications  that  detail  a  complete  model  of  the 
coherent  motions  and  their  interactions.  Initial  experiments  to  explore  ways  to  modify  the  key 
aspects  of  the  turbulence  production  mechanism  were  made.  We  attempted  to  make  change  in 
the  evolution  of  isolated  simulated  features  as  well  as  in  the  fully  turbulent  layer. 
Unfortunately  the  understanding  of  the  difficulties  expressed  by  the  statement  in  the  proposal 
"we  do  not  expect  that  optimal  controls  will  be  discovered",  were  underestimate.  The  problem 
of  quantifying  the  sensitivity  of  specific  turbulence  mechanisms  to  controls  by  making  direct, 
fixed,  local  modifications  turns  out  to  be  much  more  difficult  than  anticipated.  To  make 
progress  with  this  approach,  which  we  now  feel  will  only  work  on  simulations,  one  needs 
significantly  more  resources.  We  knew  this  going  in,  and  therefore  focused  on  the  pocket 
vortex  as  the  most  likely  mechanism  that  controls  could  be  applied  to.  Using  the  structural 
picture  of  Falco  (see  1991  Proc  Roy  Soc.,  and  earlier  references  given  below),  we  focused  on 
reducing  the  spreading  rate  of  this  vortex. 

The  two  mechanisms  that  are  operative  and  that  contribute  to  the  spreading  of  the  pocket 
vortex  are  viscous  shear,  and  inviscid  stretching  (due  to  the  presence  of  the  wall).  We 
attempted  to  change  both  the  shear  and  the  strength  of  the  'mirror  image'  effect  in  the 
experiments  performed. 

The  proposed  phased  momentum  bursts  proved  too  difficult  for  us  to  implement.  The  bottom 
line  was  that  for  all  three  'stages'  of  the  proposed  program  -  the  simulation,  the  water  tunnel 
LIPA  runs,  and  the  wind  tunnel  boundary  layer  vorticity  probe  studies  --  changes  were 
visually  observed  to  occur  (often  with  the  smallest  perturbations  we  could  generate),  but  the 
changes  that  occurred  always  seemed  to  result  in  different  evolutions  from  instant  to  instant 
even  though  the  same  perturbation  was  applied.1  This  frustrated  our  attempts  to  quantify  the 
changes.  With  the  scatter  in  the  response  so  large,  we  would  need  to  analyze  literally 
thousands  of  LIPA  grids.  Since  during  this  time  period  we  were  evaluating  grids  by  hand,  this 
was  simply  not  possible  to  do.  In  addition,  during  the  course  of  doing  this  we  discovered  from 


Iwith  hindsight  this  is  not  so  surprising  since  the  purtubations  are  applied  to  a  specific 
location,  but  the  events  to  be  modified  pass  that  loaction  with  randon  positioning  in  the 
turbulent  boundary  layer.  However,  in  the  ring/wall  simulations  we  had  reasonable  hope  to  do 
better. 


92  5  14  119 


92-12998 


2 


our  wind  tunnel  experiments  that  the  length  of  hot-wire  record  necessary  to  get  convergence 
of  vorticity  statistics  in  an  unperturbed  boundary  layer  was  so  long  that  convergence  was  poor 
in  a  natural  boundary  layer  (JFM  219, 1990).  The  1/2  million  data  points  that  we  discovered 
were  needed  for  5%  accuracy  effectively  made  the  use  of  LIPA  for  this  purpose  unrealistic. 

(It  should  be  noted  that  this  signal  length  was  100  times  longer  than  previously  reported  for 
free  shear  layers.)  Thus,  the  additional  significantly  greater  variance  observed  when 
modifications  where  introduced  would  also  make  it  essentially  impossible  to  see  the  effect  of 
the  controls  —  even  in  air.  Thus,  the  only  way  we  could  follow  this  line  of  attack  was  with  a 
drag  balance,  which  we  did  not  have. 

This  understanding  led  to  a  different  approach.  We  went  back  to  flow  visualization  in  the 
fully  turbulent  boundary  layer  in  ?;t,  and  applied  some  controls  two-dimensionally.  This  has 
three  advantages.  We  get  area  iniurmation,  we  get  Lagrangian  information  and  we  apply  the 
perturbation  in  a  similar  way  to  a  large  number  of  pocket  events.  The  controls  applied  were 
attempts  to  address  both  the  viscous  and  the  inviscid  mechanisms  that  lead  to  vortex 
stretching.  Small  separated  regions  of  the  order  of  the  viscous  sublayer  in  thickness  and  of  the 
order  of  hundreds  of  viscous  lengths  long  were  made  using  micro  rearward  facing  steps.  Flow 
visualization  of  the  formation  of  pockets  was  studied.  Results  (the  Thesis  of  S.-C.  Yoo,  see 
below)  showed  that  the  number  of  pockets  has  decreased  significantly. 

Perhaps  the  most  successful  inroads  we  have  made  into  understanding  how  to  determine  to 
sensitivity  of  the  turbulent  boundary  layer  production  mechanisms  is  summarized  by  the 
scaling  shown  in  figure  1  ( taken  from  the  proceedings  of  the  Royal  Soc  Lond  A,  336  see 
below).  I  asked  the  question,  "statistically  what  should  the  momentum  and  time  scales  of  a 
momentum  perturbation  be  if  I  am  going  to  precisely  alter  the  essential  coherent  motions?  It 
would  have  to  be  the  characteristic  momentum  and  time  scales  of  these  motions.  If  I  could 
find  those  scales,  and  if  they  are  the  important  Reynolds  stress  producing  scales,  then 
boundary  layer  statistics  such  as  the  turbulence  intensities  and  Reynolds  stresses  would  show 
similarity  (i.e.  collapse)  when  normalized  by  them.  Figure  1  shows  that  we  have  found  those 
scales  (in  the  inner  region).  This  has  been  done,  for  the  first  time,  using  coherent  motion 
ideas.  The  results  are  very  convincing.  We  now  have  a  handle  on  these  quantities.  This 
collapse  allows  us  to  predict  intensities  and  Reynolds  stress  magnitudes.  From  the  point  of 
view  of  control  we  now  know,  on  a  statistical  basis,  how  much  momentum  to  use  in  a  control 
scheme,  and  for  how  long  to  apply  it.  for  any  Reynolds  number  and  position  in  the  layer.  An 
experimental  program  could  be  initiated  again  using  these  quantities,  but  it  would  have  to  be 
applied  over  a  large  enough  area  to  get  a  statistical  result  ( basically  it  must  be  done  on  a  large  es 
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plate  attached  to  a  drag  balance).  Furthermore,  using  this  modeling  we  can  now  study  some 
aspects  of  the  sensitivity  of  the  controls  analytically.  For  example,  if  we  make  the  Typical 
Eddy  scale  smaller,  say,  by  5%,  or  the  Typical  Eddy  faster  than  nature  makes  it  in  a  normal 
boundary  layer  at  a  given  Reynolds  number,  we  can  determine  how  much  less  Reynolds  stress 
we  will  have. 


We  have  gone  a  step  further.  The  very  recent  PhD  work  of  K.  Pan  (March  1992),  which  is  an 
essential  piece  of  information,  especially  if  control  is  to  be  handled  statistically,  finally  was 
completed  in  March  1992.  It  provides  us  with  statistics  on  the  scales  and  probability  of  phase 
relationships  of  the  Typical  Eddies,  the  long  streaks  and  the  pockets  (which  give  us 
information  on  both  the  pocket  vortices  and  primary  hairpin  vortices).  If  an  area  sensing  array 
could  be  fabricated,  correctly  phased  trains  of  input  could  be  introduced,  and  statistically  we 
could  distribute  momentum  slugs  of  the  right  magnitude  and  duration  to  have  a  statistical 
match  with  the  coherent  structures.  The  results  of  tweaking  the  ’knobs"  of  this  analog 
computer  would  be  fascinating. 
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’’A  coherent  structure  model  of  the  turbulent  boundary  layer”.  Proc.  Workshop  on 
Turbulent  Wall  Layers,  ed.  R.  Panton  and  D.  Bogard,  1987. 

"What  happens  to  the  large  eddies  when  net  drag  reduction  is  achieved  by  outer  flow 
manipulators?",  (with  N.  Rashidnia).  Proc.  of  Symposium  on  Turbulent  Drag 
Reduction  by  Passive  Means,  London  1987. 

"Measurement  of  two-dimensional  fluid  dynamic  quantities  using  a  photochromic  grid 
tracing  technique",  (with  C.C.Chu),  SP1E,  Vol.  814,  Photomechanics.  1987. 

Two-dimensional  Fluid  Dynoamic  Measurements  Using  LIPA"  Proc  Image 
Processing  and  Analysis  Workshop",  ed  Brodkey  and  Guezennec,  Ohio  State  Univ. 
1988. 

"Vortex  ring/viscous  wall  layer  interaction  model  of  the  turbulence  production  process 
near  walls",  (with  C.  C.  Chu).  Exp.  in  Fluids  6  p.  305-315, 1988. 
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measurements  over  an  area",  (with  C.  C.  Chu,  M.  H.  Hetherington,  C.  P.  Gendrich). 
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"The  Turbulence  Burst  Detection  Algorithm  of  2.  Zaric’",  (with  C.  P.  Gendrich)  in 
Turbulence  Near  Walls,  ed.  S.  J.  Kline,  Hemisphere  Press,  N.Y.  1989. 
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Coherent  Structure  Dynamics  and  Suggestions  for  Improved  Devices".  AIAA  Paper 
89-1026, 1989T 

"Production  of  Turbulence  in  Boundary  Layers",  (with  J.  C.  Klewicki  and  K.  Pan)  in 
the  Proceeding  of  the  7th  Turbulent  Shear  Flows  Symposium,  Stanford  Univ., 

August,  1989. 

"Vorticity  Reid  Measurements  Using  Laser  Induced  Photochemical  Anemometry 
(LIPA)  (with  C.  P.  Gendrich  and  C.  C.  Chu).  in  the  Proceeding  of  the  7th  Turbulent 
Shear  Flows  Symposium,  Stanford  Univ.,  August,  1989. 

"Production  of  Turbulence  n  Boundary  Layers  and  Modification  of  the  Near  Wall 
Region", (With  J.  C.  Klewicki  and  K.  Pan)  in  Structure  of  Turbulence  and  Drag 
Reduction,  ed.  A.  Gyr  Springer- Verlag.  pp.59-68, 1990. 

"On  Accurately  Measuring  Statistics  Associated  With  Small  Scale  Structure  in 
Turbulent  Boundary  Layers  Using  Hot-Wire  Probes  (with  J.  Klewicki),  J.  Fluid  Mech. 
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On  the  sign  of  the  instantaneous  spanwise  vorticity  component  in  the  near  wall 
region  of  turbulent  boundary  layers,  accepted  for  publication  in  Physics  of  Fluids. 

"A  coherent  structure  model  of  the  turbulent  boundary  layer  and  its  ability  to  predict 
Reynolds  number  dependence".  Phil  Trans  Royal  Soc.  Lond.  A  336, 103-129, 1991. 

Papers  still  in  the  publication  process 

"On  the  Reynolds  Number  Dependence  of  Turbulence  Statistics  in  Boundary  Layers, 
submitted  to  Journal  of  Fluid  Mechanics. 

"Reynolds  Stress  Transport  in  Turbulent  Boundary  Layers",  to  be  submitted  to  Journal  of 
Fluid  Mechanics. 

"Organized  vortical  motions  in  the  outer  region  of  turbulent  boundary  layers",  submitted  to 
Physics  of  Fluids. 

Papers  which  have  only  had  their  ABSTRACTS  published: 

"Instantaneous  vorticity  measurements  over  an  area  in  unsteady  and  turbulent  flows" 
(with  C.  C.  Chu).  Bull.  Am.  Phy.  Soc.,  Series  II,  Vol.  32, 1987. 

"Two-point  vorticity  correlations  in  a  turbulent  boundary  layer"  (with  J.  Klewicki),  Bull. 
Am.  Phy.  Soc.,  Series  II,  Vol.  32, 1987,  pg.  2061. 

"The  bursting  process  in  turbulent  boundary  layers"  (with  C.  C.  Chu,  K.  Pan  and  J. 
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"The  role  of  outer  region  manipulator  wakes  in  turbulent  boundary  layer  drag 
reduction"  (with  N.  Rashidina),  Bull.  Am.  Phy.  Soc.,  Series  II,  Vol.  32,  pp.  2088, 1987. 

"Turbulence  Production  in  Boundary  Layers"  Proc  USC  Turbulence  Symposium,  ed. 
J.  McMichaels,  June  1988 
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"Vorticity  Field  Measurements  Using  Photoexcited  Chemicals",  Proc  USC  Turbulence 
Symposium,  e4J-  McMichaels,  June  1988. 

"Reynolds  Number  Effects  on  the  Structure  of  Turbulence  in  Boundary  Layers"  (with 
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"Eduction  of  Inner  Region  Turbulent  Boundary  Layer  Structure  for  Two-Point 
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Series  II  pg  2293. 

Need  for  High  Reynolds  Number  Experiments  Written  discussion  Symposium  of 
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series  of  lectures),  Stanford  Univ.  22-24  April  1 987. 
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Coherent  structure  model  of  turbulent  boundary  layers 


Figure  1<J-  (a)  The  rms  gtreamwise  fluctuation  (u')  distributions  over  a  wide  range  of  Reynolds 
numbers,  normalized  using  ur  and  v/ut.  eb,  x  ,  a,  Rt  —  1010,  2870  and  4850  (Klewicki  ft  Falco 
1990);  a,  =  1340(Purtell  etal.  1981);  *,R,  =  8000  (Klebanoff  1955) ;  o.a.O.R,  =  18300,28600 
and  39000  (Petrie  et  at.  1990).  (6)  The  rms  gtreamwise  fluctuation  (u')  distributions  over  a  wide 
range  of  Reynolds  numbers,  normalized  using  uT>,  and  Cv.  Symbols  are  the  same  as  in  figure  10a. 


Figure  1^  (a)  The  Reynolds  stress  <«V>  distribution  over  a  factor  of  15  in  Reynolds  number, 
normalized  using  uj,  v/ur.  ®.  x  ,  a.  Rt  —  1010,  2870  and  4850  (Klewicki  ft  Falco  1990) ;  •,  R,  —  8000 
(Klebanoff  1955);  ©,  Rt  =  14500  (Bradshaw  ft  Morrison  1986).  (6)  The  Reynolds  stress  <uV> 
distribution  over  a  factor  of  15  in  Reynolds  number,  normalized  using  uj..  and  Cr  Symbols  are  the 
same  as  in  figure  11a 

again  we  see  that  the  overlap  region  (as  defined  by  the  law  of  the  wall  for  the  mean 
flow)  is  much  better  accounted  for  by  these  length  and  velocity  scales.  In  addition, 
the  fully  turbulent  portion  of  the  boundary  layer  is  also  very  well  accounted  for.  A 
universal  power  law  relation  also  exists  for  the  Reynolds  stress,  over  approximately 
the  same  range  as  found  for  rmsu.  The  power  law  region  is  fitted  by 

The  collapse  in  the  near  wall  region,  with  about  the  same  accuracy  as  achieved 
using  wall  variables,  is  again  suggesting  that  the  te  influence  is  as  important  as  the 
wall  scales  to  within  a  few  viscous  units  from  the  wall.  The  fact  that  both  the  rmsu 
and  )  distributions,  which  are  significantly  different  in  form  (when  normalized 
by  the  classical  wall  coordinates),  can  be  collapsed  by  the  te  scales  from  the  wall 
region  through  the  fully  turbulent  region  (where  the  wall  scaling  fails)  suggests  the 
importance  of  the  te  to  boundary -layer  dynamics. 

Phil  Tram  R  Soc.  Land.  A  (1991) 


9 


